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ABSTRACT 

In a series of papers, we developed a technique for estimating the inner eccentricity 
i-C ■ in hierarchical triple systems, with the inner orbit being initially circular. However, 

for certain combinations of the masses and the orbital elements, the secular part of 



Q ■ the solution failed. In the present paper, we derive a new solution for the secular part 

\ of the inner eccentricity, which corrects the previous weakness. The derivation applies 

^ ■ to hierarchical triple systems with coplanar and initially circular orbits. The new 

formula is tested numerically by integrating the full equations of motion for systems 
with mass ratios from 10^^ to 10^. We also present more numerical results for short 
term eccentricity evolution, in order to get a better picture of the behaviour of the 
inner eccentricity. 
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INTRODUCTION 



Generally, stars have a tendency to form groups of different multiplicity, from the smallest possible (binary systems) up to 
large groups, like globular clusters with a population of the order of lO'^ stars. Observations give values for the frequency of 
^ ' multiple stars in the galactic field of up to 70%, (Kiseleva-Eggleton & Eggleton 2001 and references therein). Most of the 
, multiple systems seem to be hierarchical (e.g. Tokovinin 1997, Kiseleva-Eggleton & Eggleton 2001). 

' A hierarchical triple system consists of a binary system and a third body on a wider orbit. The motion of such a system 

can be pictured as the motion of two binaries: the binary itself (inner binary) and the binary which consists of the third body 
and the centre of mass of the inner binary (outer binary). For most hierarchical triple stars, the period ratio X is of the order 
of 100 and these systems are probably very stable dynamically. However, there are systems with much smaller period ratios, 
like the system HD 109648 with X = 22 (Jha et al. 2000), the A Tan system, with X = 8.3 (Fekel & Tomkin 1982) and the 
CH Cyg system with X = 7.0 (Hinkle et al. 1993). 

When the inner pair of a hierarchical triple stellar system is a close binary system, i.e. the separation between the 
components is comparable to the radii of the bodies, the evolution of the inner binary can depend very sensitively on the 
separation of its components and this in turn is affected by the third body. Thus, a slight change in the separation of the 
binary stars can cause drastic changes in processes such as tidal friction and dissipation, mass transfer, accretion and mass 
loss due to a stellar wind, which may result in changes in stellar structure and evolution. For instance, an eccentricity of only 
~ 0.001 can be important in a semidetached binary (Kisela, Eggleton & Mikkola 1998). A comprehensive summary of those 
topics can be found in Eggleton (2006). 

In a close binary, the orbit is expected to become circular eventually, due to tidal friction. However this is not possible when 
the binary motion is affected by the perturbations of a third companion (Mazeh & Shaham 1979, Mazeh 2008). Eggleton and 
Kiseleva (1996), based on results from numerical integrations of coplanar, prograde and initially circular orbits of hierarchical 
triple systems, derived the following empirical formula for the mean inner eccentricity: 

where A and B depend on the mass ratios. 

Equation ^ was confirmed analytically by Georgakarakos (2002). Moreover, the analytic derivation was generalized to 
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hierarchical triple stellar systems with other orbital characteristics (Georgakarakos 2003 and 2004). Subsequently, some of the 
formulae were also tested for systems with planetary mass ratios (Georgakarakos 2006). 

In the present paper, we deal with a problem that had to do with the secular part of our analytical solution and we 
manage to obtain a new formula which is free of that problem and gives improved results for certain cases. The new derivation 
applies to systems with initially circular and coplanar orbits, i.e. the type of systems we dealt with in Georgakarakos (2002). 
The new formula is tested numerically for systems with stellar and planetary mass ratios. 

Other recent work on the dynamics of hierarchical triple systems includes the work done by Krymolowski & Mazeh (1999), 
Ford, Kozinsky & Rasio (2000), Blaes, Lee & Socrates (2002) and Lee & Peale (2003). 



2 THEORY 

First, we would like to give a brief summary describing the derivation of the analytic formula (for more details, one can check 
Georgakarakos 2002). Our eccentricity expression consisted of two parts: a short period part, which varied on a time-scale 
comparable to the inner and outer orbital period, and a secular part. 

In order to obtain the short period part of the inner eccentricity, and by using Jacobian coordinates, we expanded the 
perturbing potential in terms of Legendre polynomials and we then used the definition of the Runge-Lenz vector to obtain 
a series expansion for the inner eccentricity in powers of 1/X (the period ratio X was considered to be rather large for the 
applications discussed). 

The secular part of the eccentricity was derived by using a Hamiltonian which was averaged over the short period time- 
scales by means of the Von Zeipel method. Then, by manipulating the canonical equations and by making a few approximations, 
we obtained our secular solution and eventually, we combined the short period and the secular solution in order to get a unified 
expression for the inner eccentricity (we also had to find an expression for the outer eccentric vector, in order to determine 
the required initial conditions for the secular problem). 

However, the secular solution was singular for certain combinations of the orbital elements and masses of some systems, 
and as a result, our model failed to describe the eccentricity evolution for those systems. The problem originated from the 
approximations we made about the motion of the outer pericentre, but with a more careful handling, that can now be corrected 
(see subsection 1 2. 2 |l . 

Finally, we would like to mention that an important aspect of the theory that was developed, was the combination of 
the short period and secular terms in the expressions for the eccentricities. At any moment of the evolution of the system, 
we considered that the eccentricity (inner or outer) consisted of a short period and a long period (secular) component, i.e. 
e = fishort + esec (ouc cau picture this by recalling the expansion of the disturbing function in solar system dynamics, where 
the perturbing potential is given as a sum of an infinite number of cosines of various frequencies). Thus, considering the 
eccentricity to be initially zero leads to eghort = — Cscc (initially), which implies that, although the eccentricity is initially zero, 
the short period and secular eccentricity may not be. 

We would like to point out here, that our derivation was done in the context of the gravitational non-relativistic three 
body problem with point masses. However, whenever it is necessary, it would be possible to generalize the calculation to 
include relativistic effects and the effect of treating the bodies as non-point masses, by adding the required terms (e.g. a 
post-Newtonian correction) in the equation of motion of the inner binary. 



2.1 Short period eccentricity 

First, we list the components of both the inner and outer eccentric vectors (short period). Although the derivation details can 
be found in Georgakarakos (2002), the inner eccentric vector equations in that paper contain only the two dominant terms of 
the expansion, instead of the four we give below (this is done for better accuracy purposes; a four term expansion was used 
in Georgakarakos 2003 and 2004). The expansion of the outer eccentric vector includes one more term, compared to what we 
found in Georgakarakos (2002). That extra term improves the general solution of the problem, as it provides better initial 
conditions for the secular problem. 

The components of the inner eccentric vector are: 



eii(t) = xi{t) + C, = '^^(p^2i(t) + ^PMt)+m*xip^si{t) + m,\p^S2it)]+C,,, (2) 



ei2{t) = yi(t)+Coi2 = (^Py2l{t) + j^Py22{t)+m.xiPy3l{t)+m,^Py32{t) ] +Cci2 (3) 

where 

11 9 

-Px2i(t) ~ — - cosnit -I- - cos ((3ni — 2n2)t — 2<^) + - cos ((ni — 2n2)t — 2(^) 
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1 9 

Px22(t) = -cos((3ni -2n2)t-20) + -cos((ni -2n2)t-2</>) 
6 2 

15 

-Pxsi (i) = — cos {nit + 0) 

3 45 15 

-Px32(t) = — COS ((2ni — n2)t ~ 'P) ~ cos ((2ni — Sm)^ — 30) — — cos ((4ni — 3n2)i — 30) 

11 9 
-Py2i(t) = — - sinmt + - sin ((3ni — 2n2)t — 20) — - sin ((ni — 2n2)t — 20) 

1 9 
Py22{t) = -sin((3ni - 2n2)t- 20) - -sin((ni - 2n2)i - 20) 

15 

-Py3i (t) = sin (n2t + 0) 

3 45 15 

Py32{t) = — sin ((2ni - n2)t ~ 4>) + sin ((2ni - 3n2)t ~ 30) - — sin ((4ni - 3n2)t - 30) 

and 

ru'z — rni 



{nil + rn-z) 3 Ma 

The components of the outer eccentric vector, which are obtained in a similar way to the one we used to get the components 
of the inner eccentric vector, are: 

/ \ / \ ^ TTllTfli lr3 / ,\ 1/3 . 

e2i(t) = xi{t) + Ce2i = 4 2- — tIt cos(n2t + 0) + cos((2ni - n2)t - 0) - 

(mi +m2)3M3 X3 4 A 16 

-^cos((2ni-3n2)t-30))]+Ce2i (4) 

e22(t) = j/2(t) +Ce22 = V^ll^Hl. ^_L[^sin(n2t + 0) + -^(-:^sin((2ni-n2)t-0) + 

(mi+m2)3M3X3 4 X 16 

+ ^cos((2ni-3n2)t-30))] + Ce22- (5) 

M is the total mass of the system, m and 712 are the mean motions of the inner and outer binary respectively, is the 
initial relative phase of the two binaries, i.e. the initial angle between the two Jacobi vectors r and i?, and finally, Cenfiei^^ 
Ce2i and Ce22 are constants of integration. 

We would like to mention here, that when we expanded the components of the eccentric vector as a power series in terms 
of 1/X , terms of the same order as P22 were missed out. Those terms arose from the P5 Legendre polynomial in the expansion 
of the perturbing potential and they are of the form 

105 ms 111 m* — mf 

In the best case scenario, a rough estimate gives 

P22/P51 ~ 6, 

which leads us to conclude that the impact of the P51 terms would not be significant for the parameter ranges discussed in 
this paper. This will be confirmed by numerical integrations in the subsequent sections. For that reason, and in order to avoid 
further increase of the volume of the equations, we choose not to include those terms here. 



2.2 Improved secular solution for the eccentricity 

As in our previous papers, and in order to describe the long term motion of the system, we use a Hamiltonian which is averaged 
over the inner and outer orbital periods by means of the Von Zcipel method. There are also other methods of studying the 
secular behaviour of hierarchical triple systems, e.g. see Libert & Henrard (2005), Migaszewski & Gozdziewski (2008), both 
concentrating on planetary systems. However, those methods cannot reproduce fully the results of the Von Zeipel averaging 
method (we are referring to terms that arise from the canonical transformation and which are second order in the strength of 
the perturbation). As a result, the outcome is similar for planetary mass ratios, but for systems with comparable masses, the 
extra terms that the Von Zeipel method produces can have a noticeable effect (e.g. see Krymolowski and Maaeh 1999). 
The doubly averaged Hamiltonian for coplanar orbits is (Marchal 1990): 



H = 
where 



(6) 
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Gmim2m3ag 



<3i = -o- —77; r-x(2 + 3es), (7) 

° (mi + m2)ai\,{l — e^,) 2 

15Gmim2m3(mi - m2)a|eseT , , „ 2n /qs 

Q2 = ^ — — — 1 — cos(gs -5T)(4 + 3es), (8) 



64(mi +m2)2a:*,(l 



4)' 



15 GmimzmjaleKl-ei)^ , „ 2. , „ 2 ^l 

- — ^ 3 [5(3 + 2eT) +3eT cos2(gs -5t)]. 



(mi + m-z)^ M^a^{l 



4)3 



(9) 



The subscripts S and T refer to the inner and outer long period orbit respectively, while g is used to denote longitude of 
pericentre. The first term in the Hamiltonian is the Keplerian energy of the inner binary, the second term is the Keplerian 
energy of the outer binary, while the other three terms represent the interaction between the two binaries. The Qi term 
comes from the P2 Legendre polynomial, the Q2 term comes from the P3 Legendre polynomial and the Qa term arises from 
the canonical transformation. The expression for Q3 in Marchal 1990, contains just the term which is independent of the 
pericentres. The complete expression we use here is taken from Krymolowski and Mazeh (1999). 

Using the above Hamiltonian, the secular equations of motion are the following: 

da;s 5 ex 2^,.. , g 2x . , , 2 . ^ , i^-4)^ , 25 3 + 24 
~A~ = Tr"7 — ^(l-es)'[(4 + 3es)smgT+6(a;s2/scos5T + 2/ssmpT)]-[- TTT + -5-771 ZT^si^- 

'^T^ J^D (1 — e4,)2 (1 — 3^)2 O (l — Crj,) 

-^ei)]ys + y7 ^j^^'^2-)3 [2^s cos2ffT - a;ssin25T - y(Ki + 3yi) cos 2gT + Ks(a;i + 2y|) sin 2^x1 (10) 

dt/s 5 bt 2,lr,, o 2^ . 2 M r(l-e|)i 25 3 + 2e^ ,^ 

— = "T^^: — ^(l-es)2[(4 + 3es)cossrT + 6(a;si/ssm5T + a;scosffT)] + [- ZTT + -^^T^ 

1^ J-D (1 — e^)2 (1 — e;|,)2 » ^ ^t) 

-^es)]xs + (I )3 + ys sin 2ffT - ^(t/i +3a;s)cos2gT - ys{yi + 2x1) sm2gT] (H) 

X ((a^s - J/s) cos 2srT + 2a;sys sin 2gT)] (12) 

dfiT 5 Q/3 , 2w • ^ 15 eT(l-es)5, ,2 2x . „ n .^on 

'a~ ^ TfiTi 2T2(4 + 3es)(j/scos3T -issmgx) - -^/37— — 5-(22;sys cos 2gT - (xg - yg) sm 2gT) (13) 

ar ID (^1 — e^j o (1 — e;|,)7 

where 

xs = escosgs, ys^essinps, 

1 1 ^ 

mi-m2as ^ mim2M2 ,03,1 m3 /O-s.^ j j 3 G^msa^ 

a = , (3 = 3-( — 7 = — I -( — )' and dr = 2 — 

mi+m2aT m3(mi + m2)2 M2(mi+m2)2 4 ^3 _I_ ^j^^) 2 

In Georgakarakos (2002), we assumed that remained constant, we neglected terms of order e| and and only the 
term proportional to (3 was retained in equation (|12|l . Under those assumptions, the system was reduced to one that had the 
following form: 



d£s 
dr 

dys 
dr 

dgT 
dr 



"Bys + CsinpT 

Bss — Ccos^T (14) 
A, 



where A,B, and C were constants. 

The solution of the above system was singular for A ~ B — (this is to be expected, as the above system of differential 
equations corresponds to the motion of a forced harmonic oscillator with a constant forcing frequency A and a natural 
frequency B). As a consequence, we got an overestimate of the inner eccentricity for the parameter values that satisfied 
A — B ^ 0. The source of the problem was the assumption that the outer pericentre frequency was almost constant, or in 
other words, the fact that we only kept the dominant term in the approximate expression for ^t. In fact, although the term 
proportional to a/3 in the equation for the outer pericentre is insignificant compared to the leading term in most of the cases, 
there are systems (the systems for which A~ B ~ 0) for which the a/3 term becomes comparable to the leading term, and 
hence, its inclusion in the relevant equation is necessary ( note that bt is also proportional to a/3 to leading order). 



Improved equations for HTS 5 



In order to deal with the extra term in the equation for the outer pericentre, we introduce the following two variables (as 
we did for the inner binary): 

XT ~ eTCOsgT, 3/T ~ eTsin^T- 

By using equations (12) and (13) and without making any assumption about the outer binary pericentre and the two eccen- 
tricities, we obtain expressions for xt and yx- Then, if we neglect terms of order e| and (keep in mind that we deal with 
initially circular orbits), we obtain the following system of differential equations: 



-Bys + CyT 
Bxs — Cxt 

CDys - DyT (15) 
-CDxs + DxT, 



dxs 

dr 
dys 

dr 
dxT 
"dT 
dyr 

dr 

where 

B = l + —7, C = -a and _D = /3. 
8 4 

Note that C is slightly different from the one in system (|14|) fC = (5/4)cieT in system (|14p '). 

System (|15p is a system of four first order linear ordinary differential equations. The characteristic equation of the system 

is 

a" + {D^ + + 2C^D)A^ + {DB - C^Of = (16) 
and it has the following four imaginary eigenvalues: 
Ai,2 = iifca, A3, 4 = ±ik4, 
where 



h = ^ fci4_fc2 ^ ^4 = ^\''\ ki=D'' + B^ + 2C^D and = {D + By[{D - B)^ + 4C^ D] 
The corresponding eigenvectors are 

^^^l^U,^2^,^3^,^f, i = 1,2,3,4 



with 



CC^D + \1-BD . C{B + D) 1 BC^D- {B'^ + Xl)D 



^" A, B2 + C'^D + \l ' B2 + C2D + A2 ' A, + C^D + A^ ' 

Now, since we have found the eigenvalues and the corresponding eigenvectors, we are able to obtain the general solution 
of system (|15p . which is the following: 

X = Cif7i3 + C2?723 + C3Lri4+C4C/24, (17) 

where 

X = [xs,ys,XT,yTY' , 

Uij = [Fij sin kjT, F2j cos kjT, Fsj sin kjr, cos kjr]'^ , 

U2j = [^Fij cos kjT, F2j sin kjr, —Fsj cos kjr, sin kjrf" , with j = 3, 4. 

The F coefficients are 

_CC^D-k]-BD _ C{B + D) _ 1 BC^D- {B'' - k])D 

^"^^ ^ k' B^ + C^D - fc? ' ^""^ ^ WTCPD^^' ^ fc7 B2 _^ C2JJ _ fc2 . 

and finally, the d constants, for initially circular orbits, are found to be (see second from the end paragraph of section [2}: 

Cl = -T, ^-^(yiO - y20F24), C2 = — — ^ (3:^10^34 - 2:20^14), 

-r24 — -^23 i<13-^34 — I'asru 

= — ^-^^{y2oF23 - yio), and C4 = ^r^^ — ^ „ „ (a:20-Fi3 - a;io-F'33), 

-^24 — -^23 -^13-^34 — -^33-^14 

where xio = a;i(0), yio = yi(0), 2:20 = 3:^2(0) and 1/20 = 2/2(0). 
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2.3 The new formula 

Now, if we combine the short period solution as given in subsection (2.1) and the new secular solution we derived in the 
previous section (this is done the usual way, i.e. by replacing the constants in equations ([2} and ((3| with the equations for 
xs and xt obtained from (|17[) . as the latter evolve on a much larger time-scale than the short period terms and they can 
practically be treated as constants for time-scales comparable to the inner and outer orbital period) and if we average over 
time and over the initial phase <j), we get: 



225 
^256' 



43 61_1_ 8361 2_1_ 365 1 
Y ^ y X ^ 4096 '^*3d ^ T8"X2 
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+ C'lav) 



(18) 
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3 NUMERICAL TESTING 

In order to test the validity of the formulae derived in the previous papers, we integrated the full equations of motion 
numerically, using a symplectic integrator with time transformation (Mikkola 1997). 

The code calculates the relative position and velocity vectors of the two binaries at every time step. Then, by using 
standard two body formulae, we computed the orbital elements of the two binaries. The various parameters used by the code, 
were given the following values: writing index Iwr = 1, method coefficients al = 1 and a2 = 15, correction index icor = 1. 
The average number of steps per inner binary period NS, was given the value of 60. 

For our simulations, we use the two mass ratios we defined in Georgakarakos (2006), i.e. 

Kl = and K2 = 



mi + 7712 77ll + 7712 

with 0.001 < Kl < 0.5 and 0.001 < K2 < 1000 (although our main interest is for systems with comparable masses, we also 
test the theory for some planetary mass ratios too), but we studied a few more pairs than in Georgakarakos (2006). We also 
used units such that G = 1 and mi + m2 — 1 and we always started the integrations with ai = 1. In that system of units, the 
initial conditions for the numerical integrations were as follows: 




7^1 = 0, 7^2 = 1, R 
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Figure 1. Percentage error p against initial phase 4> for Ki = 0.001 and K2 = 1000. 



3.1 Short period evolution 

First, we present some results from testing equations ^ and As in all our previous papers, we compare the averaged 
over time numerical and theoretical inner eccentricity. In addition to that, we also compare the averaged over time and initial 
phase numerical and theoretical inner eccentricity in order to get a picture of the error dependence on 0. 

Generally, the error gets larger for increasing values of K2 and, for fixed mass ratios, it goes down as the period ratio 
increases. Also, the error for a specific K2 seems to be almost independent from the variation of Ki (for a specific value of 
(j)). An interesting finding was, that for some systems, the error was seriously affected by the choice of the initial phase (j), 
especially for systems with larger K2, i.e. large ms. For instance, for K2 ~ 1000, X — 10 and (j) — 90°, the error is around 
28%; however, that error can even vanish for certain values of the initial phase. Since the formulae depend on (j>, we expected 
the error to depend on the initial phase too, but that result was a bit surprising. Figure 1 is a plot of the percentage error 
p against the initial phase for Ki — 0.001 and K2 ~ 1000. The shape and size of the curves remain almost the same as 
Ki varies. The graphs are also similar when we plot them for smaller K2- However, as K2 gets smaller, the amplitude of the 
oscillation of the error becomes smaller too, and for very small outer masses {K2 < 0.1), the curves are almost flat. 

Figures 2 show the percentage error between the averaged (top row graphs: averaged over time inner eccentricity with 
— 90°; bottom row graphs: averaged over time and initial phase inner eccentricity) numerical and theoretical inner eccen- 
tricity. The left graphs are for 0.001 < K2 < 1 while the right graphs are for 0.001 < K2 < 50. For K2 > 50, the error remains 
almost the same as for K2 ~ 50. All graphs are for Ki — 0.05. For other values of Ki, we get similar graphs, as the error does 
not vary much for the rest of the values of Ki . 

Finally, simulations that were performed including the P51 terms in equations and showed an extra improvement 
of around 1%. 

3.2 Long period evolution 

In this section, we present the results from comparing equation (|18p with the numerical solution of the full equations of 
motion. 

Each triple system was integrated numerically for = 0°— 360° with a step of 45° . After each run, e?^ was averaged over 
time using the trapezium rule and after the integrations for all finished, we averaged over by using the rectangle rule. 
The integrations were also done for smaller steps in (p (e.g. 1°), but there was not any difference in the outcome. 

These results are presented in Table 1, which gives the percentage error between the averaged numerical ef^ and equation 
(I18[l . The error is accompanied by the period of the oscillation of the eccentricity, which is the same as the integration time 
(wherever there is no period given, we integrated for an outer orbital period, i.e. 2nX, as there was not any noticeable secular 
evolution). There are four values per {KI, K2) pair, corresponding, from top to bottom, to X = 10, 20, 30, 50 respectively. 

Generally, the new formula seems to work better and to provide more consistent results than the previous one. For most 
systems that show secular evolution, our model (the averaged square eccentricity) has an error of about 10% for X = 10, 
independently of the mass ratios of the system. There are a few systems with larger error (e.g. for K2 = 0.5, but these systems 
are dominated by short period evolution) and one with a small negative value ( Ki = 0.3, K2 = 0.1 with —3% error). This 
specific system looks a bit odd in Table 1, but if there were more systems in the table, the change in the error would look 
more smooth (e.g. the error for a system with K\ — 0.2 and K2 =0.1 is 4.5%). Basically, the error we get from the short 
period solution (and which is the main source of the error in the final formula) is passed to the secular solution by means of 
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the initial conditions. As a result, our model has a longer period and some difference in the amplitude is also noticeable. As 
X increases, the period and amplitude differences go down and consequently the error goes down too. 

Regarding the new secular solution we derived in subsection (|2.2p . it appears that we have achieved our aim, i.e. to correct 
the problem that arose from the singularity A — B — in the old solution. For example, if one checks the corresponding table 
in Georgakarakos (2006), he sees that for Ki = 0.05 and K2 = 0.01, the errors for X = 10, 20, 30, 50 are -1.7, -15.9, -17, -37 
respectively. Using our new formula, the corresponding errors become 9.4, 3.9, 3, 1.6. Figures 3 demonstrate the improvement 
that is achieved by the use of our new formula. Both column graphs are plots of inner eccentricity against time. The left 
column graphs are for a system with Ki = 0.3, K2 = 0.09, X = 10 and (j) = 90°, while the parameter values for the right 
column graphs are Ki = 0.4, K2 = 0.065, X = 50 and (j) = 90°. The top graphs come from the numerical integration of the 
full equations of motion, the middle graphs are based on the old theoretical solution of Georgakarakos (2002) and the bottom 
graphs are produced by our new equations. The improvement is clearly noted (period and amplitude of oscillation). Note 
that, although X — 50 for the right column graphs, and one may expect not to have a significant error, the old equations fail 
badly. 

Finally, we would like to mention here, that , after a more careful inspection of the systems in Table 1 of Georgakarakos 
(2006), a few secular periods were corrected and some other systems that were not given a secular period, they now have one, 
as it can be seen in Table 1 of this paper. 



4 SUMMARY 

An even distant or light companion is expected to inject eccentricity to a close binary. A slight change in the binary separation 
can have drastic consequences on its subsequent evolution. Therefore, it is very useful to obtain analytical expressions about 
the amount of eccentricity injected into the binary orbit and the time-scales on which it is done. In addition, the analytical 
expressions could also be used for other purposes, such as the detection of a companion to a binary system (Mazeh and 
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Table 1. Percentage error between the averaged numerical and equation JTSj. For each {Kl,K2) pair, there are four entries, 
corresponding, from top to bottom, to X = 10, 20, 30, 50 respectively. The numbers in parentheses are the secular periods of the systems. 
In the system of units used, the inner binary period is Tin = 27r. 
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Figure 3. Inner eccentricity against time. Left column graphs: Ki = 0.3, K2 = 0.09, X = 10 and (p = 90°. Right column graphs: Ki = 0.4, 

K2 = 0.065, X = 50 and (p = 90°. The top grajjlis come from the numerical integration of the full equations of motion, the middle graphs 
show our theoretical solution based on the secular solution of Georgakarakos (2002) and the bottom graphs are based on the new secular 
solution. 



Shaham 1979), the determination of orbital elements of stellar triples or to set constrains to the masses and orbital elements 
of exoplanets. 

In a series of papers, we derived formulae for estimating the inner eccentricity in hierarchical triple systems with various 
orbital characteristics and with the inner eccentricity being initially zero. The secular part of the solution we obtained in 
those papers, failed for certain combinations of the masses and the orbital parameters of some systems. In the current paper, 
we managed to obtain a secular solution for the inner eccentricity that had no such weakness. However the derivation applies 
only to hierarchical triple systems with initially circular outer binaries. For non-circulax outer binaries, the system of the 
approximate secular differential equations is non-linear and therefore we can not obtain an improved solution as we did for 
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the circular case. The numerical integrations of the full equations of motion confirmed that our new model works well where 
the previous one failed to do so. 

Our future aim is to see whether a similar solution can be obtained for the non-coplanar formula of Georgakarakos (2004) . 
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